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ABSTRACT 

The  development  of  the  solidified  microstructure  in  metal-matrix  composites  depends  on 
complex  interactions  between  non-planar  solidification  fronts  and  multiple  particles.  The  problem 
is  multiscale  in  nature;  the  motion  of  the  particle  (under  the  action  of  a  nano-scale  disjoining 
pressure  force  and  a  micro-scale  viscous  drag  force)  is  dynamically  coupled  with  the  developing 
solidification  front  morphology,  which  is  dependent  on  a  variety  of  thermal  conditions.  Using 
computational  techniques  discussed  in  parts  I  and  II,  this  paper  seeks  to  describe  the  complicated 
nonlinear  parametric  dependencies  of  the  phenomenon.  The  effects  of  four  of  the  most  important 
parameters  in  the  particle-solidification  front  interaction  are  investigated,  i.e.  the  Hamaker 
constant,  the  particle  size,  the  thermal  conductivity  ratio  of  the  particle  to  the  melt,  and  the  solid- 
liquid  interfacial  free  energy.  By  performing  simulations  using  the  multiscale  approach  the 
dependencies  of  the  critical  velocity  on  these  four  parameters  is  clarified. 

1.  INTRODUCTION 

Particle-solidification  front  interactions  are  abundant  in  nature  and  industrial  applications. 
Examples  include:  frost  heaving  in  soils  [1],  cryobiology  [2],  electrification  of  clouds  [3],  and 
microstructural  development  in  metal-matrix  composite  manufacture  [4-24].  This  paper  is  a 
culmination  of  an  effort,  begun  in  Parts  I  and  II  [17,18]  to  develop  and  apply  a  multiscale  model 
to  study  the  interaction  of  solidification  front  with  ceramic  particles  embedded  in  the  melt,  a 
situation  that  occurs  in  the  formation  of  microstructures  in  metal-matrix  composites. 
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During  the  casting  process  of  metal  matrix  composites,  solidification  of  the  metal  proceeds  in  the 
presence  of  ceramic  particles.  A  particle  that  is  approached  by  an  advancing  solidification  front 
will  be  engulfed  if  the  solidification  velocity  exceeds  a  certain  “critical  velocity”  [7],  For 
velocities  below  the  critical  value  the  particle  is  pushed  steadily  ahead  of  the  front.  Determining 
the  conditions  under  which  a  particle  is  pushed  or  engulfed  will  aid  in  the  prediction  of  the 
overall  microstructure  of  a  system  and  enable  control  of  the  solidification  process  to  design 
desired  particle  distributions  in  metal-matrix  composite  microstructures. 

Consider  a  particle  that  is  being  approached  by  a  directionally  solidified  front  (of  a  pure  material) 
as  shown  in  Figure  1  [17,  18].  As  the  solidification  front  approaches  the  particle  distortions  in  the 
thermal  field  can  occur  when  the  particle  thermal  conductivity  {kp)  is  different  from  the 
solidifying  melt  thermal  conductivity  {ki}  [10,11].  These  distortions  will  cause  the  solidifying 
interface  to  deviate  from  planarity  and  either  form  a  convex  protuberance  (kp/ki  <  1.0)  or  a 
concave  trough  (kp/ki  >  1.0)  underneath  the  particle  [10-12].  A  protuberance  underneath  the 
particle  promotes  pushing  whereas  a  trough  underneath  the  particle  has  a  tendency  to  promote 
engulfment  the  particle  [12].  In  fact,  experiments  have  shown  that  in  cases  where  a  trough  has 
formed  {kp/ki  >  1.0)  the  particle  is  engulfed  nearly  all  of  the  time  whereas  when  a  protuberance  is 
formed  pushing  or  engulfing  can  occur  [5,6,1 1]  depending  on  the  solidification  velocity  and  other 
parameters. 

As  the  distance  between  the  solidification  front  and  the  particle  {d)  approaches  tens  of  nanometers 
a  disjoining  pressure  arises  in  the  gap.  The  disjoining  pressure  results  due  to  (repulsive) 
intermolecular  forces  that  exist  in  the  gap  between  the  two  surfaces,  namely  the  solidified 
material  and  the  ceramic  particle.  A  common  continuum  representation  for  the  disjoining 
pressure  (Ft)  is  given  by  [5]: 


n  = 


A 

67td^ 


(1) 


where  A  is  the  Hamaker  constant  and  d  is  the  distance  between  the  two  surfaces.  In  this  work,  the 
disjoining  pressure  is  defined  such  that  a  negative  Hamaker  constant  results  in  a  repulsive  force 
between  the  two  interfaces  whereas  a  positive  Hamaker  constant  results  in  an  attractive  force 
between  the  interfaces.  The  value  (and  sign)  of  the  Hamaker  constant  depends  on  the  system  in 
question,  typical  values  being  ±10’'^  Joules. 
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The  Hamaker  constant,  and  hence  the  disjoining  pressure  can  also  be  related  to  the  surface  free 
energies  of  each  of  the  phases  involved  [5]: 


A  =  -4mIqA/  (2) 

where 

^r  =  rsp-r,p-rsi  (3) 

where  ysp  is  the  interfacial  free  energy  of  the  solid(solidification  frontj-particle  interface,  yip  is  the 
liquid-particle  interfacial  free  energy,  and  ysi  is  the  solid-liquid  interfacial  free  energy.  The  Aq  is 

a  constant  with  the  dimensions  of  length  and  has  a  typical  value  of  several  molecular  diameters 
[5].  In  systems  where  Ay  is  positive,  as  the  solidification  front  gets  close  enough  to  the  particle 
the  existence  of  a  melt  layer  between  the  front  and  particle  is  energetically  favored,  in  which  a 
repulsive  disjoining  pressure  is  formed.  The  disjoining  pressure  tends  to  push  the  particle  away 
from  the  front.  The  opposite  is  true  when  the  interfacial  free  energy  difference  (Ay)  is  negative, 
i.e.  the  solidification  front  is  attracted  to  the  particle. 

The  particle  reacts  to  the  disjoining  pressure  and  begins  to  move.  Note  that  if  one  neglects  the 
volume  change  upon  solidification  when  starting  from  an  initially  quiescent  melt  there  is  no  fluid 
flow  in  the  system  until  the  particle  begins  to  move  due  to  the  disjoining  pressure.  The  motion  of 
the  rejected  particle  causes  a  drag  that  tends  to  push  the  particle  back  towards  the  solidification 
front.  Overall,  the  imbalance  between  the  drag  and  disjoining  pressure  forces  determines  whether 
the  particle  will  be  pushed  or  engulfed.  The  following  four  effects  impact  on  the  magnitudes  of 
the  drag  and  repulsive  disjoining  pressure  forces: 

1 .  Thermal  conductivity  ratio:  The  magnitude  of  the  forces  acting  on  the  particle,  particularly  the 
drag,  depends  on  the  curvature  of  the  solidification  front.  For  example,  it  has  been  shown  that  the 
drag  acting  on  the  particle  significantly  increases  as  the  kp/ki  increases  [12].  This  is  due  to  the 
large  viscous  losses  in  the  gap  as  the  solidification  front  becomes  less  convex  as  kp/ki  increases. 
This  effect  is  illustrated  in  pressure  contours  for  various  kp/ki  values  shown  in  Figure  2. 

2.  Hamaker  constant/wettability  of  the  particle:  Due  to  the  existence  of  a  disjoining  pressure  in 
the  melt  gap  between  the  particle  and  the  front  a  change  in  the  melting  temperature  of  the 
interacting  interfaces  will  result  through  the  Clausius-Clapeyron  relation  [3].  If  Ay  >  0,  a  layer  of 
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liquid  (premelted  fdm)  will  exist  between  the  solidification  front  and  the  particle  even  when  the 
temperature  is  lower  than  that  of  the  bulk  melting  temperature  of  the  solidifying  medium.  The 
solid-liquid  interface  temperature  (T,)  in  the  presence  of  such  a  premelting  film  is  given  by  [13]: 


T.=T- 

I  m 


A 

\dj 


where  X  is  an  interaction  length  scale  and  Tm  is  the  bulk  melting  temperature. 


(4) 


The  existence  of  a  premelted  film  and  the  interfacial  temperature  depression  can  have  a 
significant  effect  on  front-particle  interaction.  The  premelting  effect  alters  the  shape  of  the 
solidification  front  and  hence  changes  the  way  the  forces  behave.  If  one  assumes  premelting,  the 
solidification  front  deforms  to  accommodate  a  liquid  layer  -  hence  affecting  the  forces  that  act  on 
the  particle  (see  Figure  3)  [14].  As  shown  in  Figure  3,  the  solidifying  interface  begins  to  bend 
around  the  particle  in  order  to  accommodate  the  liquid  layer.  This  change  in  curvature  of  the  front 
leads  to  significant  changes  in  the  forces  that  act  on  the  particle. 


3.  Solid-liquid  interfacial  tension  effect  through  the  Gibbs-Thomson  condition:  The  solid/liquid 
interfacial  free  energy  of  the  solidification  front  also  affects  the  shape  and  curvature  of  the 
solidifying  interface.  The  solid/liquid  interfacial  free  energy  enters  the  problem  through  the 
modification  of  the  melting  temperature  in  Eq.  (4)  as  follows  [15,16]: 


T..=T- 


T 

m  TT  ^  fn 

pHsi 


(5) 


The  last  term  on  the  right  accounts  for  the  Gibbs-Thomson  effect.  The  y*/  is  the  solid-liquid 
interfacial  free  energy  of  the  solidification  front,  k  is  the  curvature  of  the  solidification  front,  p  is 
the  density  of  the  melt,  and  H^i  is  the  latent  heat  of  fusion  per  unit  mass.  The  Gibbs-Thomson 
effect  reduces  the  magnitude  of  the  curvature  of  the  solidification  front  (i.e.  it  tends  to  flatten  the 
front  out).  Hence,  the  forces  will  be  affected  and  the  critical  velocity  will  change. 
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4.  Particle  size  (R^):  The  particle  radius  also  plays  a  role  in  the  determination  of  the  critical 
velocity.  As  the  particle  radius  increases  the  drag  on  the  particle  increases  at  a  faster  rate  than  the 
repulsive  disjoining  pressure  force  acting  on  the  particle.  This  is  simply  due  to  the  fact  that  the 
repulsive  force  only  acts  over  a  small  portion  of  the  particle  surface  (i.e.  the  portion  where  the 
gap  is  small  enough  so  that  the  disjoining  pressure  effects  become  important),  whereas  the  fluid 
drag  acts  on  the  entire  particle  surface.  Models  and  experiments  have  both  shown  that  the  critical 
velocity  is  inversely  proportionally  to  the  particle  radius  [4-6]. 

All  of  the  work  on  particle-front  interactions  up  to  now  has  relied  on  focusing  on  one  parameter 
or  another,  mostly  in  a  steady-state  pushing  mode  [6,7,10,13,15,16].  There  have  been  a  few 
dynamic  models  [4,  8,  22]  of  particle-solidification  front  interactions  but  they  too  have  relied 
heavily  on  various  assumptions  to  make  the  problem  tractable  (e.g.  no  premelting,  no  surface 
tension  effects,  simplified  solidification  models  and  simplified  force  models).  This  study 
includes  all  of  the  main  parameters  involved  (the  particle  radius  (Rp),  the  kp/ki,  the  Hamaker 
constant,  and  the  solid-liquid  interfacial  energy)  and  illustrates  how  these  effects,  when  acting  in 
concert,  makes  for  a  highly  sensitive  dynamical  problem.  The  methodology  for  modeling  and 
solving  the  resulting  multiscale  problem  was  detailed  in  Part  I  [17].  In  Part  II  [18]  the  mechanism 
for  determination  of  the  critical  velocity  was  elucidated.  Here,  a  detailed  parameteric  study  is 
performed  to  clarify  the  dependence  of  the  critical  velocity  on  each  of  the  above-listed  four  main 
parameters  that  influence  the  particle-front  interactions.  In  addition,  the  physics  behind  the 
dependency  of  the  critical  velocity  on  the  parameters  listed  above  is  also  clarified. 

2.  COMPUTATIONAL  SETUP 

The  computational  setup  is  shown  in  Figure  4.  The  solidification  front  moves  upward  and 
approaches  the  particle.  In  order  to  capture  all  of  the  essential  physics  of  the  problem  in  a 
multiscale  setting,  the  computational  domain  needs  to  be  divided  into  two  different  regions.  In 
the  “outer”  region  (comprising  of  the  particle  and  the  flow  and  thermal  field  surrounding  it)  the 
full  2D  flow  and  heat  transport  equations  are  solved.  In  the  “inner”  region  (comprising  of  the  gap 
between  the  particle  and  the  front)  a  lubrication  equation  is  solved  by  making  a  thin  gap 
assumption  {d  «  Rp).  The  solutions  of  the  “inner”  and  “outer”  regions  are  then  coupled  through 
a  “matching”  region  by  supplying  boundary  conditions  to  each  region.  The  fully  coupled  solution 
in  the  two  regions  proceeds  in  tandem  over  each  time  step.  The  forces  on  the  particle  are 
evaluated  using  the  composite  solution  (i.e.  inner  and  outer  region  solutions)  and  the  particle  is 
moved  in  Lagrangian  fashion.  The  calculation  of  the  dynamics  of  the  particle-front  interaction 
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takes  into  account  all  of  the  most  important  physical  mechanisms  in  a  multiscale  setting.  The 
details  of  the  numerics  are  discussed  in  parts  I  and  II  [17,18]  and  will  not  be  discussed  in  this 
paper. 

3.  RESULTS  AND  DISCUSSION 

3.1.  Physical  Description  of  the  Pushing-Engulfing  Transition  Based  on  Front  Curvature 
Arguments 

The  curvature  of  the  solidification  front  directly  beneath  the  particle  (henceforth  referred  to  as  the 
“tip”  of  the  solidification  front)  has  a  significant  impact  on  the  critical  velocity.  Figure  5  shows 
the  evolution  of  the  particle-front  system  for  a  two  different  cases,  one  where  the  particle  is 
approached  by  an  initially  planar  front  (Figure  5(a),  (c))  and  is  pushed  ahead  of  it  and  the  other 
where  the  particle  becomes  engulfed  by  the  front  (Figure  5(b),(d)). 

Case  1  (Pushing):  Figures  5  (a)  and  (c)  respectively  show  the  particle  velocity  (v^)  versus  time  (t) 
and  the  ratio  of  the  tip  curvature  to  (also  referred  to  as  the  dimensionless  tip  curvature  /c,*) 
versus  time  for  a  case  where  the  particle  is  pushed  by  the  front,  i.e.  Jip  =  1  micron,  ^=-8x10“'^ 
Joules,  /l=3.5xl0''°  meters,  y^/  =0,  and  a  solidification  velocity  (v^)  of  200  p/s.  As  shown  in 
Figure  5  (c),  for  this  case  the  tip  curvature  starts  out  at  a  value  approximately  half  the  particle 
radius  and  then  settles  to  a  steady  mean  value  not  long  after  the  start  of  the  particle-front 
interaction.  Correspondingly,  as  can  be  seen  from  Figure  5  (a),  the  particle  velocity  significantly 
increases  and  begins  to  oscillate  at  a  value  commensurate  with  the  solidification  front  velocity, 
i.e.  it  enters  a  pushing  mode. 

Case  2  (Engulfment):  Figures  5  (b)  and  (d)  denote  the  particle  velocity  versus  time  and  the  tip 
curvature  versus  time  for  the  same  system  but  for  a  somewhat  higher  solidification  velocity  at 
which  the  particle  is  engulfed  by  the  front,  i.e.  Vj  =  215  p/s.  In  contrast  with  Case  1,  in  Figure  5 
(d),  the  curvature  of  the  tip  of  the  solidification  front  also  starts  out  around  2.0  but  then  reaches  a 
maximum  and  begins  to  decrease.  The  tip  curvature  continues  to  decrease  until  it  reaches  a  value 
of  around  -1.  This  indicates  that  the  shape  of  solidification  front  is  conforming  to  the  shape  of  the 
particle. 
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The  decrease  in  the  tip  curvature  in  Figure  5(d)  is  due  to  the  premelting  effect.  Since  there  is  a 
tendency,  according  to  the  premelting  theory,  for  a  liquid  layer  to  persist  in  the  gap  between  the 
particle  and  the  front,  the  tip  of  the  solidification  front  begins  to  flatten  and  assume  a  concave 
shape  in  conformity  with  the  shape  of  the  particle.  The  flattened  solidification  front  tip  results  in 
a  larger  drag  since  there  will  be  larger  viscous  losses  in  the  melt  flowing  into  the  narrow  gap. 
The  large  drag  causes  the  particle  to  slow  down,  as  can  be  seen  in  Figure  5  (b).  With  further 
flattening,  the  curvature  of  the  tip  becomes  such  that  the  drag  outweighs  the  disjoining  pressure 
forces.  Once  the  drag  becomes  the  dominant  force,  the  particle  begins  to  slow  down  to  an 
eventual  stop  leading  to  engulfment.  On  the  other  hand,  in  the  case  where  the  particle  is  pushed 
ahead  of  the  front,  the  particle  achieves  a  high  enough  velocity,  relative  to  the  tip  velocity,  so  that 
the  gap  thickness  never  reaches  a  point  where  premelting  effects  dominate.  Then  the  interface  tip 
curvature  does  not  decrease  significantly  (as  seen  in  Figure  5(c))  and  hence  the  drag  force  does 
not  dominate  the  repulsive  disjoining  force 


A  further  aspect  of  the  effect  of  solidification  front  curvature  needs  to  be  considered,  i.e.  the 
curvature  of  the  front  prior  to  contact  with  the  particle.  Figure  6  plots  the  tip  curvature  versus 
time  for  both  the  pushing  and  engulfment  cases  that  were  described  above.  It  is  seen  that  the  tip 
curvature  in  the  pushing  (i.e.  Vj  =  200  p/s.  Figures  5(a),  (c))  case  starts  out  as  a  higher  value  than 
in  the  engulfment  case  (v^  =  215  p/s.  Figures  5(b),  (d)).  In  other  words,  the  pushing  case  has  a 
more  pronounced  protuberance  underneath  the  particle  than  the  engulfment  case  even  before  the 


particle  begins  to  interact  with  the  front. 


For  the  higher  Stefan  number,  Ste  = 


k,G, 

- (where 

PH,V, 


Gi  is  the  thermal  gradient  in  the  liquid),  corresponding  to  the  lower  solidification  velocity  the 
front  senses  the  insulating  particle  earlier  and  deviates  from  planarity  before  approaching  the 
particle.  For  the  higher  front  velocity,  since  the  Stefan  number  is  low,  the  front  remains  nearly 
planar  as  it  approaches  the  particle. 


In  determining  whether  the  particle  is  pushed  or  engulfed  by  the  front,  there  is  a  delicate  balance 
between  the  development  of  a  premelted  layer,  the  speed  of  the  solidification  front,  and  other 
parameters  (particle  size,  thermal  conductivity,  solid-liquid  interfacial  energy)  that  affect  how  the 
solidification  front  evolves  in  time.  Each  of  these  effects  impact  on  the  curvature  of  the  front 
directly  beneath  the  particle,  which  as  discussed  above  correlate  with  whether  the  particle  is 
pushed  or  engulfed  by  the  front.  Therefore,  a  parametric  study  is  performed  using  the  numerical 
methodology  developed  in  Parts  I  and  II  to  study  how  the  critical  velocity  depends  on  the  most 
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significant  solidification  process  and  material  parameters.  In  order  to  obtain  the  critical  velocity 
for  each  parameter  value  several  runs  were  made  to  identify  points  in  the  parameter  space  where 
the  system  shows  a  transition  from  pushing  to  engulfment  and  the  boundary  between  the  two 
behaviors  is  demarcated.  To  condense  data,  only  those  points  that  fall  on  the  boundary  between 
pushing  and  transition  are  shown  in  the  figures.  In  general,  the  number  of  cases  to  be  run  was 
decreased  considerably  by  judicious  choice  of  parameter  values  and  by  observing  the  trends  seen 
for  each  parameter  value.  The  final  results  only  are  presented  in  the  following  parametric  study  in 
order  to  make  a  concise  presentation  of  the  data  and  insights. 

3.2.  Critical  velocity  as  a  function  of  particle  size  (Rp ) 

First,  the  dependence  of  the  critical  velocity  on  the  particle  radius  (Rp)  is  studied.  Figure  7  shows 
the  solidification  velocity  far  from  the  front  (v^),  which  in  directional  solidification  terminology  is 
Vpuii,  (the  velocity  with  which  the  sample  is  “pulled”),  versus  Rp.  Other  parameters  are  kp/ki  = 
0.01,  A=  -8x10'^^  Joules,  1=3.5x10''°  meters,  and  y*/  =  0.  Several  cases  were  run  for  each  Rp  and 
the  points  plotted  demarcate  the  transition  from  pushing  to  engulfment.  The  solid  line  represents 
a  curve-fit  to  the  transition  between  pushing  and  engulfment  behavior  that  describes  the  variation 
of  the  critical  velocity  with  the  particle  radius.  From  Figure  7  it  is  observed  that  there  is  an 
inverse  relationship  between  the  critical  velocity  and  the  particle  radius  with  the  curve-fit 
providing  the  expression:  =201 1 (where  Vcr  is  in  p/s  and  Rp  is  in  microns).  This 

computed  dependency  agrees  with  the  experimentally  observed  relationship  Vcr  Rp‘  [4,5]. 

3.3.  Critical  velocity  as  a  function  of  particle  thermal  conductivity  (kp/ki ) 

The  general  effect  of  kp/ki  on  the  critical  velocity  is  well  established  in  the  particle-solidification 
front  interaction  literature.  It  is  well  known  and  it  can  be  physically  explained  that  when  kp/ki  <  1 
particle  pushing  can  result,  while  for  kp/ki  >  1  engulfment  is  the  most  likely  result.  However,  it  is 
difficult  to  precisely  quantify  the  specific  effect  of  the  kp/ki  experimentally  since  changing  the 
kp/ki  would  require  changing  the  materials  under  study  and  consequently  the  Hamaker  constant,  1. 
From  theoretical  work,  Shangguan  et  al.  [6]  and  Catalina  et  al.  [4]  suggest  that  the  critical 
velocity  is  inversely  proportional  to  the  kp/ki.  More  specifically,  they  obtain  a  critical  velocity  of 


(6) 


V, 


CR 


la^A/ 

/ 

64*3juRj, 


where  is  the  molecular  diameter  and  is  on  the  order  of  Angstroms  in  size  and  jU  is  the  dynamic 
viscosity  of  the  melt. 


Note  that  in  determining  the  critical  velocity  obtained  in  Eq.  (6)  simplifying  assumptions  were 
made  to  render  the  analysis  tractable.  In  that  work  the  premelting  effect  was  not  included  in  Eq. 
(6)  which  required  an  ad  hoc  cut-off  gap  thickness  to  be  invoked  in  order  to  distinguish  between 
pushing  and  engulfment  (see  references  [4,6]  for  discussion  of  a  gap  thickness  cut-off).  The 
present  modeling  and  solution  approach  allow  for  the  removal  of  limitations  inherent  in  the 
analysis  presented  in  [6].  Therefore,  the  above  expression  for  the  dependency  of  critical  velocity 
on  thermal  conductivity  ratio  is  re-examined  with  the  current  multiscale  approach. 

Figure  8  compares  the  variation  of  the  critical  velocity  with  respect  to  the  kp/ki  (with  A  =  -8x10’’^ 
Joules,  1  =3.5x10''°  meters,  Rp  =  \  micron,  and  y^;=0)  using  the  current  method  to  that  of  the 
previously  suggested  inverse  relationship  due  to  Shangguan  et  al.  [6].  From  Figure  8,  in 
agreement  with  Shangguan  et  al.  [6],  the  critical  velocity  decreases  as  the  kp/ki  increases.  This 
dependency  can  be  explained  from  the  physics  of  the  problem,  in  that  larger  kp/ki  causes  the 
solidification  front  to  flatten  out  with  the  front  eventually  becoming  concave  underneath  the 
particle  when  kp/ki  becomes  greater  than  one  (see  Figure  2).  The  drag  therefore  increases  as  kp/ki 
increases,  which  tends  to  promote  engulfment.  Since  engulfment  is  more  likely  to  occur  when 
kp/ki  increases  the  critical  velocity  decreases  with  increasing  kp/ki.  However,  Figure  8  shows  that 
the  functional  form  of  the  relationship  between  the  critical  velocity  and  the  kp/ki  is  quite  different 
from  that  given  by  Shangguan  et  al.  [6]  in  Eq.  (6).  Note  that  the  functional  relationship  indicated 
by  the  present  calculations  takes  into  account  the  delicate  interplay  between  the  drag  force  and 
the  disjoining  pressure  force.  The  strength  of  these  forces  is  highly  dependent  on  the  shape  of  the 
solidification  front.  The  kp/ki  greatly  influences  the  shape  of  the  solid-liquid  interface.  However, 
the  kp/ki  value  is  not  the  only  factor  that  influences  the  shape  (in  particular  the  curvature)  of  the 
solidification  front.  Premelting  and  Gibbs-Thomson  effects  also  play  a  role.  In  the  work  of 
Shangguan  et  al.  [6],  premelting  effects  were  not  included  and  hence  the  solid-liquid  interface 
shape  was  influenced  solely  by  kp/ki  effects. 
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3.4.  Critical  velocity  as  a  function  of  particle  wettability  (Hamaker  constant  (4)) 


As  mentioned  before,  the  Hamaker  constant  is  related  to  the  interfacial  free  energy  difference 
between  the  three  phases  (solid,  liquid,  and  particle)  through  Eq.  (2).  As  pointed  out  for  kp/ki  in 
Section  3.3,  experimentally  isolating  the  effect  of  the  Hamaker  constant  on  the  critical  velocity 
also  proves  to  be  difficult.  However,  a  comparison  can  still  be  made  with  theoretical  models, 
particularly  those  of  Shangguan  et  al.  [6]. 

Figure  9  plots  the  solidification  velocity  far  from  the  front  (v^)  versus  the  Hamaker  constant  (A) 
for  systems  where  Rp  =  1  micron,  =  0,  and  kp/k/  =  0.01.  It  is  seen  that  as  A  becomes  more 
negative  the  critical  velocity  increases.  This  is  to  be  expected,  since  as  the  value  of  the  disjoining 
pressure  in  Eq.  (1)  becomes  more  negative,  the  repulsive  force  between  the  particle  and  the 
solidification  front  becomes  greater  which  tends  to  promote  pushing.  From  Eq.  (1),  it  would 
appear  likely  that  the  Hamaker  constant  affects  the  critical  velocity  in  a  linear  fashion,  since  the 
disjoining  pressure  is  linearly  related  to  the  Hamaker  constant.  This  is  also  suggested  by  the 
model  presented  by  Shangguan  et  al.  [6].  In  fact  the  present  calculations  deviate  only  slightly 
from  linearity.  The  reason  for  the  slight  deviation  is  due  to  the  fact  that  changing  the  Hamaker 
constant  also  changes  the  premelting  parameter,  X,  through  the  equation  [13]: 


;i  = 


A 


(7) 


where,  p  is  the  density  of  the  melt  and  Hsi  is  the  latent  heat  of  fusion  of  the  melt  per  unit  mass. 


The  disjoining  pressure  in  the  premelted  layer  and  the  interface  temperature  are  related  via  the 
Clausius-Claperyon  equation.  The  change  in  the  solidification  front  morphology  due  to  the 
depression  of  interface  temperature  causes  the  forces  acting  on  the  particle  to  be  altered  and  this 
effect  is  reflected  in  the  deviation  from  linearity  of  the  curve  in  Figure  9.  In  the  work  of 
Shangguan  et  al.  [6],  no  premelting  was  assumed  and  therefore  a  linear  dependence  was  found. 
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3.5.  Critical  velocity  as  a  function  of  solid-liquid  surface  tension  (Vsi ) 


In  the  cases  computed  so  far  the  interfacial  temperature  on  the  melting  temperature  was  not 
included.  In  reality,  the  curvature  of  the  solidification  front  affects  the  temperature  at  which  the 
solidification  front  solidifies  through  the  Gibbs-Thomson  effect,  viz.: 


T,  =  T  -  T  - 

^  i  ^  m  tt  ^ 

pH,i 


A 

\dj 


(8) 


where  is  the  bulk  melting  temperature  and  /r  is  the  solidification  front  curvature. 


Non-dimensionalizing  by  GiRp,  d  by  dc,  and  /rby  gives: 


t:  =  -T/c  - 


A 


dd* 


G,Rp 


(9) 


Where  F  = 


/  sl^m 

GpRlpH, 


is  a  non-dimensional  capillarity  parameter  for  the  solid-liquid  interface. 


Interfacial  tension  at  the  solid-liquid  interface  hinders  the  formation  of  a  bump  underneath  the 
particle.  Figure  10  illustrates  how  the  critical  velocity  varies  as  a  function  of  F.  In  these  cases 
kp/ki  =  0.01,  A  =  -8x10'^^  Joules,  1  =  3.5x10’'°  meters,  and  =  1.0x10’°  meters.  From  Figure  10, 
it  is  seen  that  as  F  increases,  the  critical  velocity  decreases.  The  reason  for  this  is  that  the 
convexity  of  the  interface  underneath  the  solidification  front  is  less  pronounced  as  F  increases, 
i.e.  higher  curvatures  are  discouraged  by  capillary  effects.  Since  the  solidification  front  is  flatter 
in  the  case  of  higher  F,  the  pressure  losses  in  the  gap  are  higher  and  hence  the  drag  is  higher. 
Since  the  drag  force  tends  to  promote  engulfment,  for  higher  F  the  particle  will  more  likely  be 
engulfed,  and  hence  the  critical  velocity  is  smaller.  Figure  10  also  indicates  that  the  critical 
velocity  is  nearly  a  linear  function  of  F. 
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4.  CONCLUSIONS 


The  curvature  of  the  tip  of  solidification  front,  and  its  evolution  in  time,  has  a  direct  effect  on 
whether  the  particle  is  pushed  or  engulfed.  A  high  positive  (convex)  curvature  tends  to  promote 
pushing  since  the  fluid  drag  decreases  as  the  curvature  increases.  The  curvature  of  the  interface, 
however,  is  a  complex  time-dependent  function  of  the  thermal-fluid  transport  combined  with 
interfacial  phenomena  such  as  premelting  layers  and  the  resulting  pressure  and  temperature  in  the 
layers,  the  surface  energies  of  the  interacting  interfaces  etc.  This  computational  study  has 
clarified,  by  means  of  a  detailed  parameteric  analysis  the  following  effects  of  the  most  important 
parameters  in  particle-front  interactions: 

1)  Effects  of  pull  speed  in  the  directional  solidification  process:  As  v,  increases  the 
difference  between  the  diffusion  rate  and  the  solidification  rate  decreases  (i.e.  the  Stefan 
number  decreases)  and  hence  the  distance  at  which  the  particle  is  felt  by  the  front  is 
decreased.  This  causes  the  solidification  front  to  remain  flat  until  the  presence  of  the 
particle  begins  to  distort  the  temperature  field.  This  is,  arguably,  a  small  effect  unless  v, 
is  very  high.  The  higher  the  Vj  the  smaller  the  curvature  of  the  front  as  it  approaches  the 
particle.  This  fact,  along  with  the  fact  that  the  particle  velocity  may  not  have  time  to 
catch  up  with  and  outstrip  the  high  solidification  velocity  makes  the  particle  more  prone 
to  engulfment. 

2)  Effects  of  the  particle  velocity  (vn)  relative  to  the  tip  velocity  (v,):  If  the  tip  velocity  is 
high  enough  that  the  particle  is  not  pushed  away  sufficiently  from  the  front,  the  effects  of 
premelting  in  the  thin  melt  gap  take  over;  then  the  curvature  of  the  tip  will  decrease  and 
the  particle  will  begin  to  experience  a  large  drag.  This  drag  further  impedes  the  motion 
of  the  particle  until  eventual  engulfment  ensues. 

3)  Effects  of  particle  thermal  conductivity  (k/ki):  When  kp/ki  is  different  from  unity  the 
particle  distorts  the  thermal  field  of  the  system  which  causes  a  deformation  of  the 
solidification  front.  The  smaller  the  kp/ki,  the  larger  the  critical  velocity. 

4)  Effect  of  the  Hamaker  constant:  The  Hamaker  constant.  A,  affects  the  solidification  tip 
curvature  through  its  effect  on  the  premelted  layer.  If  a  layer  of  liquid  is  energetically 
favored  to  exist  between  the  particle  and  the  solidification  front,  the  curvature  of  the  tip 
will  change  to  accommodate  such  a  layer.  The  critical  velocity  increases  as  A  becomes 
more  negative  due  to  the  larger  repulsive  force  that  acts  on  the  particle. 
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5)  Solid-liquid  interfacial  energy  Effects  {F)  :  The  interfacial  free  energy  of  the  solid-liquid 
interface,  y^i  (or  /’in  non-dimensional  terms),  controls  the  curvature  of  the  solidification 
front.  Increasing  F  tends  to  enhance  the  smoothness  of  the  solidification  front.  Since  the 
interface  flattens,  the  critical  velocity  decreases  as  /’increases. 

The  above  effects  working  in  a  coupled  and  dynamic  fashion  make  the  problem  of  particle- 
solidification  front  interactions  a  very  challenging  one  to  analyze  from  a  theoretical  standpoint. 
The  computational  methodology  developed  in  Part  I  [17]  has  enabled  the  detailed  parametric 
study  presented  in  this  paper.  The  results  obtained  provide  insights  and  quantification  of  the 
dependencies  of  the  critical  velocity  on  the  most  important  thermo-physical  properties  in  the 
particle-front  interaction  system. 
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Figure  1.  (a)  Schematic  of  the  overall  system,  (b)  Schematic  of  the  zoomed-in  box  in  (a) 
illustrating  the  different  mechanisms  involved  in  the  particle-solidification  front  interaction.  As 
the  solidification  front  approaches  the  particle,  the  disjoining  pressure  starts  to  push  the  particle 
(if  Ay  >  0).  Fluid  then  flows  in  to  replenish  the  gap  which  results  in  a  drag  force. 
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a)  (b) 


c) 


c) 

Figure  2:  Interface  shape  and  pressure  contours  for  cases  where  the  solidification  front 
approaches  a  particle  with  the  properties:  a)  kp/ki  =  0.01,  b)  kp/ki  =  1.0,  and  c)  kp/ki  =  30.0.  The 
indicated  pressure  is  non-dimensional. 
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a)  (b) 


Figure  3:  Premelting  effect  on  the  solidification  front  (interfaces  are  shown  by  bold  curves  and 
the  contours  are  of  temperature  are  also  shown  in  the  figures):  (a)  Overall  view  and  (b)  Zoomed- 
in  vies  of  the  interface.  The  premelting  effect  leads  to  a  loss  of  planarity  of  the  interface. 
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Figure  4.  Initial  system  setup  for  the  computations  performed. 
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a)  b) 


c)  d) 


Figure  5:  Relationship  between  interface  curvature,  velocity  and  the  pushing/engulfment  of  a 
particle  when  approached  by  a  solidification  front,  (a)  Case  1  -  Particle  velocity  versus  time  for 
the  pushing  case  kp/ki  =  0.01,  Rp  =  1.0x10'*’  m,  J=-8xl0'^^  Joules,  A.=3.5xl0''°  meters,  v,  =  200 
microns/second,  (b)  Case  2-  particle  velocity  versus  time  for  the  case  where  the  particle  is 
engulfed  by  the  front;  kp/ki  =  0.01,  Rp  =  l.OxlO'*’  m,  J=-8xl0''^  Joules,  }^=3.5xl0''°  meters,  Vj  = 
215  microns/second,  (c)  and  (d)  Non-dimensional  curvature  versus  time  corresponding  to  (a)  and 
(b)  respectively. 
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Figure  6:  Non-dimensional  curvature  versus  time.  For  the  solid  line  (engulfed  case):  kp/ki  = 
0.01,  Rp  =  l.OxlO'*’  m,  J=-8xl0''^  Joules,  }t=3.5xl0''°  meters.  Vs  =  215  microns/second.  For  the 
dotted  line  (pushed  case):  kp/ki  =  0.01,  Rp  =  1.0x10''’  m,  J=-8xl0'^^  Joules,  }t=3.5xl0''®  meters, 

=  200  microns/second. 
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Rp  (microns) 

Figure  7:  Identification  of  the  critical  velocity  dependency  on  the  particle  radius  (Vj  vs.  The 
other  parameters  are  held  fixed  at  kp/ki  =  0.01,  ^  =  -8.0x10''^  J,  1  =  3.5x10'*°  m,  and  jsi  =  0.0. 


Figure  8:  Dependence  of  the  critical  velocity  on  the  ratio  of  thermal  conductivity  of  the  particle 
to  the  melt,  v,  vs.  kpiki .  The  other  parameters  are  held  fixed  at  Rp=  \  micron,  A  =  -8.0x10’'^  J,  X 
=  3.5x10’'°  m,  and  y,/  =  0.0. 
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Figure  9:  Dependence  of  the  critical  velocity  on  the  Hamaker  constant,  Vj  vs.  A.  The  other 
parameters  are  held  fixed  at  k/ki  =  0.01,  Rp  =  1.0x10''^  m,  and  ysi  =  0.0. 


Figure  10:  Dependence  of  the  critical  velocity  on  the  solid-liquid  interface  tension,  vs.  F.  The 
values  of  the  other  parameters  are  held  fixed  at  kp/ki  =  0.01,  Rp  =  1.0x10'^  m,  ^=-8x10’’^  Joules, 
J^=3.5xl0'^°  meters. 


